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ABSTRACT 


The  goal  of  this  proposal  is  to  develop  a  new  NDT  concept  and  theory  based  on  Lamb  waves,  in  which  damage  can 
be  detected  “without  relying  on  past  reference  data.”  In  particular,  a  NDT  methodology  will  be  formulated  for 
detecting  defects  such  as  crack,  corrosion  and  delamination  commonly  found  in  aircraft.  To  the  Pi’s  best  knowledge, 
no  other  group  or  literature  has  addressed  NDT  issues  from  this  new  viewpoint.  The  proposed  concept  of  “reference- 
free”  damage  diagnosis  attempts  to  address  some  of  the  fundamental  technical  hurdles  that  have  to  be  overcome 
before  NDT  systems  can  be  deployed  for  continuous  monitoring  of  aircraft.  The  intent  of  the  present  approach  is  not 
to  disregard  useful  prior  knowledge  or  information  gained  but  to  relax  explicit  dependency  on  baseline  data  so  that 
the  proposed  methodology  can  be  made  more  attractive  for  aircraft  monitoring  during  its  normal  operation.  By 
incorporating  the  proposed  approach  to  the  conventional  NDT  techniques,  it  is  envisioned  that  more  reliable  damage 
diagnosis  can  be  achieved. 

First,  the  theoretical  framework  for  the  proposed  reference-free  diagnosis  will  be  developed  based  on  the  polarization 
characteristics  of  lead  zirconate  titanate  (PZT)  wafer  transducers.  Then,  features  that  are  sensitive  to  damage  but 
insensitive  to  ambient  variations  of  aircraft  in  the  field  will  be  extracted  by  comparing  two  instantaneous  collected 
signals  without  referencing  to  past  baseline  data.  Here,  active  sensing  devices  such  as  lead  zirconate  titanate  will  be 
used  to  generate  and  measure  guided  waves  in  metallic  structures.  Once  damage-sensitive  features  are  identified, 
statistical  classifiers  will  be  developed  to  establish  decision  boundaries  without  using  prior  reference  data  and  to 
minimize  false  indications  of  damage.  Finally,  realistic  environmental  and  operational  conditions  that  in-service 
aircraft  is  subject  to  will  be  explicitly  considered. 

This  strategy  offers  a  potential  for  a  significant  breakthrough  in  SHM/NDT  practice  via  the  integration  of  active 
sensing,  smart  materials,  statistical  pattern  recognition  techniques  and  new  theoretical  development  for  damage 
diagnosis  that  has  not  been  attempted  to  date.  To  the  primary  investigator’s  best  knowledge,  no  other  group  or 
literature  has  addressed  NDT  from  this  point  of  view:  “baseline-free”  damage  diagnosis.  This  field  deployable  NDT 
technique  will  be  unique  because  (1)  damage  diagnosis  can  be  accomplished  without  relying  on  prior  reference  data, 
(3)  the  monitoring  system  will  have  its  own  “intelligence”  that  autonomously  converts  real-time  sensor  data  to 
damage  diagnosis  information  without  requiring  data  interpretation  by  the  end  users,  (4)  the  proposed  NDT  paradigm 
can  be  applied  to  complex  geometries  such  as  welded/bolted  connections  where  defects  often  occur,  and  (5) 
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environmental  and  operational  variations  that  in-service  structures  experience  will  not  affect  damage  diagnosis, 
minimizing  false  indications  of  damage. 


1.  Introduction 


The  goal  of  this  proposal  is  to  develop  a  new  NDT  concept  and  theory  based  on  Lamb  waves,  in  which  damage  can 
be  detected  “without  relying  on  past  reference  data.”  In  particular,  a  NDT  methodology  will  be  formulated  for 
detecting  defects  such  as  crack,  corrosion  and  delamination  commonly  found  in  aircraft.  To  the  Pi’s  best  knowledge, 
no  other  group  or  literature  has  addressed  NDT  issues  from  this  new  viewpoint.  The  proposed  concept  of  “reference- 
free”  damage  diagnosis  attempts  to  address  some  of  the  fundamental  technical  hurdles  that  have  to  be  overcome 
before  NDT  systems  can  be  deployed  for  continuous  monitoring  of  aircraft.  The  intent  of  the  present  approach  is  not 
to  disregard  useful  prior  knowledge  or  information  gained  but  to  relax  explicit  dependency  on  baseline  data  so  that 
the  proposed  methodology  can  be  made  more  attractive  for  aircraft  monitoring  during  its  normal  operation.  By 
incorporating  the  proposed  approach  to  the  conventional  NDT  techniques,  it  is  envisioned  that  more  reliable  damage 
diagnosis  can  be  achieved.  Specific  objectives  of  the  project  are: 

■  Task  1:  Understanding  Time  Reversibility  and  PZT  Polarization  Characteristics:  This  study  will  explore  the 
PZT  polarization  characteristics  and  time  reversal  acoustics  for  online  monitoring  of  defects  in  aircraft  structures. 
In  this  task,  theoretical,  numerical  and  experimental  investigations  will  be  conducted  to  better  understand  the  time 
reversibility  of  Lamb  waves  and  the  polarization  characteristics  of  lead  zirconate  titanate  (PZT)  wafer  transducers. 

■  Task  2:  Extraction  of  Reference-Free  Features:  Based  on  the  advanced  understanding  of  time  reversibility  of 
Lamb  waves  and  PZT  polarization  characteristics,  this  task  will  explore  the  PZT  polarization  characteristics  and 
time  reversal  process  to  extract  baseline-free  features  that  are  sensitive  to  defects  but  insensitive  to  operational  and 
environmental  variations  of  aircraft. 

■  Task  3:  Instantaneous  Damage  Diagnosis:  Once  damage-sensitive  features  are  identified,  a  theoretical 
framework  for  statistical  damage  diagnosis  will  be  developed  to  fully  automate  the  decision-making  procedure. 
Special  attention  will  be  paid  to  instantaneously  determine  decision  boundaries  for  damage  diagnosis  without 
relying  on  previously  established  thresholds  and  to  minimize  false  alarms  of  damage. 

■  Task  4:  Environmental  and  Operational  Variations:  Realistic  environmental  and  operational  conditions,  such  as 
temperature  changes  and  ambient  loading,  that  in-service  aircraft  is  subject  to,  will  be  explicitly  taken  into 
consideration  through  laboratory  specimen  tests  and  flight  tests  of  aircraft.  In  addition,  practical  instrumentation 
issues  such  as  variations  in  PZT  transducer  size,  bonding  condition  and  alignment  will  be  investigated. 

■  Task  5:  Extension  to  More  Complex  Structures:  Once  the  basic  concept  of  the  reference-free  NDT  technique  is 
established,  this  concept  will  be  extended  for  more  complex  components  found  in  aircraft.  For  instance,  this  task 
will  investigate  the  effects  of  stringers,  holes,  and  bolts  in  composite  lap  joints,  on  the  proposed  approach. 


2.  Theoretical  Background 

The  development  of  the  new  NDT  technique  proposed  in  this  project  was  based  on  the  use  of  PZT  wafer  transducers 
and  Lamb  waves  propagating  along  thin  plate  structures.  In  this  part  of  the  report,  the  effect  of  PZT  poling  direction  on 
Lamb  wave  propagation  is  discussed  first.  Then,  the  damage  detection  procedure  using  instantaneously  measured 
Lamb  wave  signals  is  introduced  briefly. 

2.1  The  effect  of  PZT  poling  directionality  on  Lamb  wave  propagation 

In  this  section,  it  is  investigated  how  the  phase  of  a  Lamb  wave  mode  changes  depending  on  (I)  the  poling  directions 
of  exciting  and  sensing  PZT  wafer  transducers  and  (2)  whether  a  wafer  transducer  is  attached  either  on  the  top  or 
bottom  surface  of  a  plate.  For  illustration,  it  is  assumed  four  identical  PZT  wafer  transducers,  labeled  as  “A”,  “B”, 
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“C”,  and  “D”,  are  attached  to  a  plate  as  shown  in  Figure  1  (a).  The  arrows  indicate  positive  poling  directions  of 
PZT  transducers.  PZTs  A  and  D  are  placed  exactly  at  the  same  position  but  on  the  other  side  of  the  plate.  PZTs  B  and 
C  are  positioned  in  a  similar  fashion.  Furthermore,  it  is  assumed  that  a  narrowband  toneburst  signal  is  applied  as  an 
input,  and  the  driving  frequency  is  chosen  such  that  only  the  fundamental  symmetric  (Sq)  and  anti-symmetric  (Aq) 
modes  are  generated.  In  this  paper,  the  term  of  “positive  bending”  is  used  when  the  positively  polarized  side  of  the 
PZT  is  subjected  to  tensile  strain.  On  the  other  hand,  the  PZT  is  subjected  to  negative  bending  when  the 
negatively  polarized  side  of  the  PZT  is  subjected  to  tensile  strain.  The  positive  bending  produces  a  “positive” 
output  voltage  while  the  negative  bending  results  in  a  “negative”  output  voltage  value. 

When  PZT  A  is  excited,  the  Sq  and  Aq  modes  are  generated  and  measured  at  PZTs  B  and  C  [1].  In  an  ideal 
condition,  the  amplitude  and  arrival  time  of  the  Sq  modes  measured  at  PZTs  B  and  C  should  be  identical.  In  addition, 
both  PZTs  B  and  C  would  be  subjected  to  positive  bending  because  of  the  symmetric  nature  of  the  Sq  mode. 
Because  both  PZTs  B  and  C  are  subject  to  the  positive  bending,  the  phase  as  well  as  the  amplitude  and  arrival  time  of 
the  So  mode  measured  at  these  PZTs  is  identical.  As  far  as  the  Aq  mode  is  concerned,  PZT  C  is  subjected  to  the 
negative  bending  although  PZT  B  still  undergoes  the  positive  bending.  Therefore,  the  Aq  modes  measured  at  PZTs  B 
and  C  are  out-of-phase.  However,  when  the  poling  direction  of  PZT  C  is  switched  [Figure  2  (a)],  PZTs  B  and  C  will 
produce  out-of-phase  Sq  modes  and  in-phase  Aq  modes  [Figure  2  (b)  and  Figure  3  (b)]. 

In  the  following  section,  the  effect  of  PZT  poling  directionality  is  further  investigated  so  that  the  mode 
conversion  due  to  crack  formation  can  be  extracted  from  the  measured  Lamb  wave  signals. 


(a)  Test  configuration  I  with  collocated  PZTs  with  the 
opposite  poling  directions 


(b)  The  So  mode  produces  the  same  bending  for 
PZTs  B  and  C  while  the  Aq  mode  results  in  the 
opposite  bending. 


Figure  1.  The  effect  of  the  PZT  poling  directions  on  the  phases  of  the  Sq  and  Aq  modes  (Configuration  I) 


(a)  Test  configuration  II  with  all  PZTs  with  the  same 
poling  directions 


(b)  The  So  mode  produces  the  opposite  bending  for 
PZTs  B  and  C  while  the  Ao  mode  results  in  the 
same  bending. 


Figure  2.  The  effect  of  the  PZT  poling  directions  on  the  phases  of  the  So  and  Ao  modes  (Configuration  II) 
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(a)  So  and  Aq  modes  measured  from  configuration  I  (b)  So  and  Ao  modes  measured  from  configuration 
in  Figure  1  (a):  So  modes  in-phase  &  Aq  modes  H  in  Figure  2  (a):  So  modes  out-of-phase  &  Aq 
out-of-phase  modes  in-phase 

Figure  3.  A  schematic  comparison  of  the  So  and  Aq  modes  measured  from  Configurations  I  and  II  shown  in 
Figure  I  (a)  and  Figure  2  (a),  respectively:  AB  (a  dash  line)  and  AC  (a  solid  line)  denote  the  response 
signals  measured  at  PZTs  B  and  C  when  a  toneburst  input  is  applied  at  PZT  A. 


2.2.  Detection  of  crack  induced  mode  conversion  using  a  PZT  poling  direction 

In  this  section,  the  PZT  polarization  characteristic  is  further  advanced  so  that  the  mode  conversion  due  to 
crack  formation  can  be  detected  without  using  any  prior  baseline  data.  First,  the  effect  of  a  crack  on  Lamb 
wave  modes  is  described.  If  Lamb  waves  propagating  along  a  thin  plate  with  a  uniform  thickness 
encounter  a  discontinuity  such  as  a  sudden  thickness  variation  of  the  plate,  some  portion  of  the  waves  are 
reflected  at  the  discontinuity  point  and  others  are  transmitted  through  it.  When  a  So  mode  arrives  at  the 
discontinuity  as  shown  in  Figure  4,  the  transmitted  wave  is  separated  into  So  and  Ao  modes  (denoted  as 
So/So  and  Ao/So,  respectively).  In  a  similar  manner,  an  Aomode  is  also  divided  into  So  and  Ao  modes  (So/Ao, 
and  Ao/Ao).  Similarly,  the  reflected  waves  are  split  to  So  and  Ao  modes.  This  phenomenon  is  called  mode 
conversion  [2]. 

In  order  to  fully  investigate  the  sign  changes  of  Lamb  wave  modes  due  to  mode  conversion,  sign 
notations  are  defined  first.  In  Figure  5,  sign  notations  for  the  So  and  Ao  modes  are  defined  schematically. 
As  shown  in  Figure  5  (a),  the  So  and  Ao  modes  are  defined  to  be  positive  when  they  cause  the  deformed 
shape  of  the  specimen’s  top  surface  to  be  convex.  On  the  other  hand,  the  So  and  Ao  modes  are  called 
negative  when  the  deformed  shape  of  the  top  surface  becomes  concave  [Figure  5  (b)].  In  Figure  6,  it  is 
shown  how  the  signs  of  Lamb  wave  modes  are  determined  as  they  transmit  through  a  crack.  From 
numerical  simulations  and  experiments,  it  has  been  shown  that  the  sign  of  a  single  (non-converted)  mode 
is  not  affected  by  crack  formation.  That  is,  a  positive  So  mode  always  produces  a  positive  So/So  mode, 
and  a  negative  Ao  mode  generates  a  negative  Ao/Ao  mode,  respectively. 

On  the  other  hand,  the  signs  of  newly  generated  modes  (Ao/So  and  So/Ao)  can  be  altered 
depending  on  the  characteristics  of  a  discontinuity  that  the  launching  Lamb  mode  is  passing 
through.  Although  the  signs  of  these  converted  modes  cannot  be  determined  without  knowing  the  detailed 
characteristics  of  the  discontinuity,  certain  relationships  among  these  converted  modes  can  be  revealed. 
For  instance,  if  a  positive  So  mode  creates  a  positive  Ao/So  mode,  a  positive  Ao  mode  also  produces  a 
positive  So/Ao  mode.  That  is,  the  signs  of  the  Ao/So  and  So/Ao  modes  should  be  always  identical. 
This  is  based  on  the  reciprocity  of  signals  AB  and  BA  [3].  Here,  signal  AB  denotes  the  response 
signal  measured  at  PZT  B  when  the  excitation  is  applied  at  PZT  A.  In  order  for  signals  AB  and  BA  to  be 
identical,  the  shape,  amplitude  and  phase  of  the  Ao/So  mode  in  signal  AB  should  be  identical  to  those  of  the 
So/Ao  mode  in  signal  BA.  In  addition,  the  sign  of  the  Ao/So  mode  created  from  a  positive  So  mode 
should  be  always  opposite  to  that  of  the  Ao/So  mode  generated  from  a  negative  So  mode. 

In  Figure  7  (a)  and  Figure  7  (b),  the  phases  of  the  So  and  Ao  modes  in  signals  AB  and  CD  are  compared 
when  the  specimen  is  in  an  intact  condition.  As  for  signal  AB,  when  PZT  A  generates  positive  So  and  Ao 
modes,  both  modes  produce  positive  bending  in  PZT  B  as  shown  in  the  upper  part  of  Figure  7(a).  On  the 
other  hand,  PZT  C  in  signal  CD  creates  a  negative  So  mode  and  a  positive  Ao  mode  because  the 
poling  direction  of  PZT  C  is  opposite  to  PZT  A  as  shown  in  the  lower  part  of  Figure  7  (a). 
However,  because  the  poling  directions  of  PZT  B  and  D  are  also  opposite  each  other,  the  negative 
So  mode  will  produce  positive  bending  in  PZT  B  and  the  positive  Ao  mode  produces  positive  bending  as 
well.  As  a  consequence,  the  phases  of  the  So  and  Ao  modes  measured  at  PZT  D  are  always  identical  to 
those  measured  at  PZT  B.  Therefore,  when  the  plate  is  in  a  pristine  condition  and  four  identical  PZTs 
are  instrumented  as  shown  in  Figure  7  (a),  it  is  concluded  that  signal  AB  becomes  identical  to  signal 
CD  as  shown  in  Figure  7  (b). 
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However,  signal  AB  is  no  longer  identical  to  signal  CD  when  there  is  a  crack  between  PZTs  A  and  B 
(or  PZTs  C  and  D)  as  shown  in  Figure  7  (c)  and  Figure  7  (d).  As  for  signal  AB,  the  So/Aomode  arrives  at 
PZT  B  earlier  than  the  Aq/Sq  mode  when  the  notch  is  located  closer  to  PZT  A  than  PZT  B  (assuming  that 
the  So  mode  travels  faster  than  the  Aq  mode).  Therefore,  the  So  mode  is  followed  by  the  So/Ao,  Ao/So, 
and  Ao  modes  in  signal  AB.  Based  on  the  sign  convention  in  Figure  6,  both  the  So/Ao  and  Ao/So  modes  are 
positive  so  that  the  signs  of  all  modes  are  the  same  in  signal  AB.  On  the  other  hand,  the  So  mode  is 
followed  by  Ao/So,  So/Ao,  and  Ao  modes  in  signal  CD  because  the  So/Ao  mode  arrives  at  PZT  D  later 
than  the  Aq/Sq  mode.  In  this  case,  a  negative  Ao/So  mode  and  a  positive  So/Aomode  are  created  because  of 
a  negative  So  mode  and  a  positive  Aomode,  respectively  [Figure  6].  Note  that  signal  CD  has  two  So  modes, 
So/Ao  and  So,  and  their  signs  are  opposite.  Also,  two  Ao  modes  in  signal  CD,  Ao/So  and  Ao,  have  opposite 
signs  as  well. 

In  Figure  7  (d),  signals  AB  and  CD  are  drawn  considering  not  only  the  arrival  times  of  each  modes  but 
also  their  relative  signs  (or  phases).  Note  that,  while  the  So  and  Ao  modes  in  Figure  7  (d)  are  in-phase,  the 
So/Ao  and  Aq/Sq  modes  in  signals  AB  and  CD  are  fully  out-of-phase.  Therefore,  the  additional  modes 
generated  by  a  notch  can  be  extracted  simply  by  subtracting  signal  AB  from  signal  CD  as  shown  in  Figure 
7  (d).  In  addition  to  signals  AB  and  CD,  it  can  be  readily  shown  that  the  So  and  Ao  modes  in  signals  AC 
and  BD  are  identical  while  the  So/Ao  and  Ao/So  modes  are  fully  out-of-phase.  Because  this  approach 
relies  only  on  the  comparison  of  two  instantaneous  signals  obtained  at  the  current  state  of  the 
system  rather  than  comparison  with  previously  recorded  reference  data,  it  is  expected  that  this 
approach  may  reduce  false  alarms  of  defects  due  to  changing  operational  and  environmental  conditions  of 
the  system. 
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Figure  4.  A  schematic  diagram  of  mode  conversion  and  reflection  due  to  a  discontinuity  on  a  plate 
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Figure  5.  Sign  definitions  of  the  So  and  Aq  modes  traveling  in  a  plate  specimen 
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Figure  6.  Relative  phase  information  of  the  So  and  Aq  modes  during  the  mode  conversion  process  due  to  crack 
formation 
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(a)  An  intact  plate  with  the  PZT  configuration  II  shown  (b)  Comparison  of  signals  AB  and  CD  without  a 

notch:  the  So  &  Ao  modes  are  identical 


in  Figure  2  (a) 
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(c)  A  damaged  plate  with  the  PZT  configuration  II  (d)  Comparison  of  signals  AB  and  CD  with  a 
shown  in  Figure  2  (a)  notch:  the  So  &  Ao  modes  are  identical,  but  the 

S0/A0&  Ao/So  modes  are  out-of-phase 

Figure  7.  Extraction  of  the  additional  Lamb  wave  modes  generated  by  a  notch  by  comparing  signals  AB  and 
CD  (Ao/So  mode  denotes  an  Ao  mode  converted  from  a  So  mode  when  it  passes  through  a  crack.  So/Ao  is 
defined  similarly). 

3.  Project  Activities 

3.1.  Extraction  of  reference-free  features 


In  the  beginning  of  the  project,  the  idea  of  extracting  converted  Lamb  wave  modes  due  to  damage  using  strategically 
placed  PZT  transducers  was  validated  first.  For  the  verification,  not  only  numerical  simulation  but  also  experiment 
was  conducted. 


3.1.1.  Numerical  simulation 


Using  COMSOL  3.3a  Multiphysics  software  (www.comsol.com).  Lamb  wave  propagation  in  a  two  dimensional 
aluminum  plate  was  simulated  using  the  combination  of  plain  strain,  piezo  plain  strain,  and  electrostatics  modules  in 
COMSOL  software.  The  length  of  the  plate  was  70cm,  and  its  thickness  was  6  mm.  Four  identical  PZTs  with  10 
mm  length  and  0.508  thickness  were  attached  to  the  plate  model  as  shown  in  Figure  8.  Note  that  PZTs  A  and  D  were 
collocated  but  on  the  other  side  of  the  plate  with  the  same  poling  direction.  PZTs  B  and  C  were  placed  in  a  similar 
fashion.  The  parameter  values  used  in  the  numerical  simulation  are  listed  in  Table  1. 

Figure  9  (a)  illustrates  that  signals  AB  and  CD  were  almost  identical  and  this  well  matches  with  the 
theory.  Once  a  notch  of  3  mm  depth  and  1  mm  width  was  introduced  100  mm  away  from  PZT  A  toward  PZT  B, 
signal  AB  became  different  from  signal  CD  as  a  result  of  the  mode  conversion  induced  by  the  crack  as 
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illustrated  in  Figure  9  (b).  In  the  following  experimental  results,  the  outcome  of  this  numerical  simulation  is  further 
substantiated. 


^100  mm|^ 

0.508  mm  71  mm  width  ^ 


6  mm 


Poling  direction 


— zuu  riirn  ► 

300  mm- 
700  mm- 


-200  mm— ► 


Figure  8.  Dimension  of  an  aluminum  plate  used  in  numerical  simulation 


Table  1.  Parameters  used  in  numerical  simulation 


Parameters 

Values 

Exciting  frequency 

150  kHz 

a  (Mass  damping  coefficient) 

10'^ 

P  (Stiffness  damping  coefficient) 

0 

Sampling  rate 

4  MS/s 

Relative  tolerance 

10'"^ 

Absolute  tolerance 

10-10 

Maximum  BDF  order 

2 

Mesh  size  (Mapped  mesh) 

1  mm  X  1  mm  maximum 

(a)  Signals  AB  and  CD  without  a  notch  (b)  Signals  AB  and  CD  with  a  notch 

Figure  9.  Comparison  of  signals  AB  and  CD  (or  signals  AC  and  BD)  when  there  is  no  notch  Comparison  of  signals 
AB  and  CD  (or  signals  AC  and  BD)  with  a  notch  of  3  mm  depth  and  1  mm  width 


3.1.2.  Experimental  results 

To  further  examine  the  proposed  reference-free  NDT  technique,  experimental  tests  have  been  conducted  on  an 
aluminum  plate.  The  overall  test  configuration  and  the  test  specimen  are  shown  in  Figure  10.  The  data  acquisition 
system  was  composed  of  an  arbitrary  waveform  generator  (AWG),  a  high-speed  signal  digitizer  (DIG),  a  low 
noise  preamplifier  (LNP)  and  a  multiplexer.  The  dimension  of  the  plate  was  122  cm  x  122  cm  x  0.6  cm,  and 
four  PSI-5A4E  type  PZT  wafer  transducers  (1.0  cm  x  1.0  cm  x  0.0508  cm)  were  mounted  in  the  middle  of  the 
plate.  PZTs  A  and  D  were  collocated  and  attached  on  the  other  side  of  the  plate,  and  PZTs  B  and  C  were  mounted  in  a 
similar  fashion.  The  PZTs  were  attached  so  that  their  poling  directions  were  identical  to  the  configuration 
shown  in  Figure  2  (a).  PZTs  A  and  B  (or  PZTs  C  and  D)  were  0.52  m  apart  each  other.  In  this  experiment,  the  PZT 
transducers  were  attached  to  either  the  top  or  the  bottom  surface  of  the  plate  with  commercial  cyanoacrylate  adhesive. 

Using  the  14-bit  AWG,  a  toneburst  signal  with  a  10  peak-to-peak  voltage  and  a  driving  frequency  of  150 
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kHz  was  generated  and  applied.  First,  PZT  A  in  Figure  10  (b)  was  excited  by  this  input  waveform.  Then,  PZT  A 
generated  elastic  waves  and  the  response  was  measured  at  PZT  B.  When  the  waves  arrived  at  PZT  B,  the  voltage 
output  from  PZT  B  was  amplified  by  the  LNP  with  a  gain  of  twenty  and  measured  by  the  DIG.  The  sampling  rate  and 
resolution  of  the  DIG  were  20  MS/sec  and  16  bits,  respectively.  In  order  to  improve  the  signal-to-noise  ratio,  the 
forwarding  signals  were  measured  twenty  times  and  averaged  in  the  time  domain.  After  the  forwarding  signals  from 
PZT  A  to  PZT  B  (signal  AB)  were  measured,  the  same  process  was  repeated  by  measuring  signals  AC,  BD  and  CD. 
The  entire  experimental  process  without  averaging  took  less  than  10  seconds. 


(a)  A  data  acquisition  system 


D  =  0.52  m 


T.=  6  mm 


PZT  A 


PZTB 


L=  1.22  m 


(b)  A  square  aluminum  plate 


,  (PZT  bih  the  back)  (PZT  C  in  the  back) 

is>  V 


Figure  10.  Testing  configuration  for  detecting  a  crack  on  an  aluminum  plate 


In  Figure  II,  signals  AB  and  CD  is  instantaneously  measured  from  the  pristine  condition  of  the  specimen  and  shown. 
Figure  11  (b)  shows  an  initial  difference  between  signals  AB  and  CD  even  in  the  absence  of  a  crack  due  to 
imperfections  in  PZT  alignment,  size  and  bonding  condition.  Next,  a  3  mm  (depth)  x  1  mm  (width)  x  60  mm  (length) 
notch  was  introduced  between  PZTs  A  and  B  (or  PZTs  C  and  D).  The  notch  was  located  150mm  away  from 
PZT  A  toward  PZT  B.  As  a  consequence,  two  additional  modes  due  to  mode  conversion  appeared  between  the 
existing  Sq  and  Aq  modes  as  shown  in  Figure  12  (a).  Comparison  of  Figure  11(b)  and  Figure  12(b)  clearly 
demonstrates  the  appearance  of  the  additional  modes  due  to  crack  formation. 

However,  it  is  not  clear  in  practice  how  large  the  signal  difference  should  be  before  it  leads  to  warning  of 
damage.  To  tackle  this  issue,  a  new  instantaneous  damage  diagnosis  technique  is  proposed  in  section  3.2. 
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(a)  Signals  AB  and  CD  without  notch 
Figure  11.  Comparison  of  signals  AB  and  CD  without  notch 


(b)  Difference  between  signals  AB  and  CD 
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(a)  Signals  AB  and  CD  with  a  3  mm  notch 


0.06  O.OS  0.1  0.12  0.14  0.16  0.18  0.2 

Time  (ms) 

(b)  Difference  between  signals  AB  and  CD 


Figure  12.  Comparison  of  signals  AB  and  CD  with  a  3  mm  depth  notch 


3.2.  Instantaneous  damage  diagnosis 

In  the  previous  section,  it  has  been  shown  that  damage-sensitive  features  can  be  identified  from 
instantaneously  measured  Lamb  wave  signals.  However,  a  statistical  damage  diagnosis  is  still  required  to 
differentiate  signals  due  to  damage  from  error  signals  produced  by  variations  in  PZT  size,  alignment,  and 
bonding  condition.  In  this  section,  attention  has  been  paid  to  determining  decision  boundaries  for  damage 
diagnosis  without  relying  on  previously  established  thresholds. 

3.2.1.  Damage  classification  using  instantaneously  measured  Lamb  wave  signals 

Signals  AB  and  CD  are  indistinguishable  when  there  is  no  crack  as  shown  in  Figure  6.  This  is  based  on  the 
assumption  that  all  PZT  transducers  are  identical  and  PZTs  A  and  D  or  (PZTs  B  and  C)  are  perfectly  collocated.  In 
practice,  these  assumptions  cannot  be  fully  satisfied  because  of  variations  in  PZT  size,  alignment  and 
bonding  condition  [4].  This  imperfection  in  PZTs  may  generate  initial  differences  in  signals  AB  and  CD  even  in  the 
absence  of  crack  and  lead  to  positive  false  alarms. 

Here,  a  damage  detection  scheme  is  developed  based  on  the  premise  that  mode  conversion  produces 
signal  differences  between  signals  AB  and  CD  that  are  bigger  than  the  initial  differences  due  to  PZT  imperfection. 
The  proposed  technique  takes  advantage  of  not  only  signals  AB  and  CD  but  also  signals  AC  and  BD  to  extract  mode 
conversion  in  the  presence  of  variations  of  the  PZT  size,  alignment  and  bonding.  The  uniqueness  of  the  proposed 
damage  classification  scheme  is  that  threshold  values  for  damage  classification  is  obtained  using  only  signals 
instantaneously  measured  from  the  current  state  of  the  system  without  relying  on  predetermined  threshold  values. 

In  Figure  13,  signal  AB  is  schematically  drawn  based  on  the  sign  definition  presented  in  section  2.2 
assuming  a  notch  is  closer  to  PZT  A.  From  Figure  13,  it  can  be  shown  that  signal  AB  is  a  simple  superposition  of 
signals  So,  MCl,  MC2  and  Aq.  Signal  So  indicates  a  time  signal  that  contains  only  the  So  mode  but  its  length  is 
identical  to  that  of  signal  AB.  Signal  MCl,  signal  MC2  and  signal  Ao  are  defined  in  a  similar  fashion.  MCl  and  MC2 
represent  the  first  and  second  arrivals  of  Lamb  wave  modes  created  by  mode  conversion,  respectively.  Note  that  MCl 
and  MC2  could  be  either  So/Ao  or  Ao/So  modes  depending  on  the  relative  position  of  crack  and  the  actuating 
and  sensing  PZTs  used  for  the  signal  measurement.  For  instance,  MCl  denotes  a  So/Ao  mode  in  signal  AB  when  a 
notch  is  closer  to  PZT  A  than  PZT  B.  This  is  because  the  So/Ao  mode  arrives  at  PZT  B  earlier  than  the  Ao/So  mode. 
Similarly,  signal  MC2  in  signal  AB  represents  the  Ao/So  mode. 

Similar  to  signal  AB,  signals  AC,  BD,  and  CD  can  be  also  expressed  as  combinations  of  individual 
Lamb  wave  signals  [Figure  14].  In  Figure  14,  these  four  signals  are  schematically  shown  emphasizing  the  relative 
phases  of  individual  Lamb  modes  among  these  signals.  For  instance,  the  Sq  and  MCl  modes  in  signal  AC  are  out-of- 
phase  compared  to  these  modes  in  signal  AB.  Therefore,  signal  AC  can  be  obtained  by  flipping  signals  So  and  MCl 
and  summing  up  signals  Sq,  MCl,  MC2  and  Ao  all  together.  Signals  BD  and  CD  are  related  to  signals  So, 
MCl,  MC2  and  Ao  in  similar  manners.  Based  on  these  observations,  the  relationship  between  the  signals  that  can 
be  measured  (signals  AB,  AC,  BD  and  CD)  and  the  individual  Lamb  signals  (signals  So,  MCl,  MC2  and  Ao)  can  be 
obtained  as  follows: 
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Because  of  measurement  errors  and  variations  in  PZT  size,  alignment  and  bonding  condition,  error  terms  are  added  in 
Equation  (1).  Here,  e^  is  an  error  signal  in  the  measured  signal  AB  that  is  superimposed  to  the  exact  signal 
AB,  and  e^^.,  Cgj^,  and  are  defined  similarly.  Cg^^  represent  error  signals  in  each  decomposed 

Lamb  mode  signals.  While  Equation  (1)  illustrates  how  different  combinations  of  individual  Lamb  mode  signals 
constitute  each  measurable  signal.  Equation  (2)  below  shows  how  each  individual  mode  signal  can  be  extracted 
from  the  measured  signals  AB,  AC,  BD  and  CD.  For  instance,  signal  So  can  be  extracted  by  adding  signals  AB  and 
CD  and  subtracting  signals  AC  and  BD.  Similarly,  signal  MCi  can  be  isolated  by  adding  signals  AB  and  BD  and 
subtracting  signals  AC  and  CD. 
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So  far,  it  has  been  assumed  that  a  notch  is  located  closer  to  PZT  A  than  PZT  B  if  it  exists.  In  reality,  the  location 
of  the  notch  is  unknown  in  advance.  When  the  notch  is  formed  closer  to  PZT  B,  the  MCl  and  MC2  modes  in  signals 
MCl  and  MC2  are  switched.  That  is,  the  MC2  mode  appears  in  the  estimated  signal  MCl  instead  of  the  estimated 
signal  MC2.  (Note  that  the  definition  of  the  MC2  mode  is  the  second  arrival  of  the  converted  mode,  and  it  appears  in 
the  Aq  region  of  signal  MCp)  From  these  observations,  it  is  found  that  Equation  (2)  can  be  still  used  regardless  of  the 
location  of  the  notch.  In  fact,  it  can  be  used  to  decide  whether  the  notch  is  closer  to  PZT  A  or  B.  For  example,  if  a 
mode  found  in  signal  MCl  is  located  in  the  So  region,  the  notch  is  located  closer  to  PZT  A.  Conversely,  if  the  MC2 
mode  is  found  in  the  Aq  region  of  signal  MCl  instead  of  the  So  region,  it  confirms  that  the  notch  is  closer  to  PZT  B. 
Therefore,  once  the  existence  of  mode  conversion  is  determined,  the  location  of  the  defect  can  be  also  determined  by 
checking  the  location  of  the  mode  in  signal  MCp 

According  to  Equation  (2),  the  measured  signals  can  be  separated  into  individual  Lamb  mode  signals  So,  MCl, 
MC2  and  Ao.  Note  that  only  signals  MCl  and  MC2  contain  Lamb  modes  created  by  damage  (MCl  and  MC2, 
respectively).  Ideally,  signals  MCl  and  MC2  should  be  zeros  in  the  absence  of  damage.  In  practice,  however,  due  to 
Cmci  and  eMC2  that  are  superimposed  to  the  exact  signals  MCl  and  MC2,  signals  MCl  and  MC2  may  not  be  zeros 
even  without  damage.  Since  Cmci  and  eMC2  cannot  be  extracted  from  signals  MCl  and  MC2  using  Equation  (2)  only, 
it  is  challenging  to  determine  whether  additional  modes  in  these  signals  are  due  to  mode  conversion  or  PZT 
imperfections.  To  tackle  this  issue,  a  new  damage  classifier  is  developed  to  determine  if  the  non-  zero  responses  in 
signals  MCl  and  MC2  are  due  to  mode  conversion  or  simply  due  to  initial  errors. 

The  damage  identification  scheme  proposed  in  this  study  starts  from  the  decomposition  of  individual 
Lamb  mode  signals  from  measured  signals  using  Equation  (2).  Next,  the  maximum  absolute  amplitude  of  the 
error  components,  Cso  and  Cao,  in  signals  So  and  Aq  are  estimated.  Because  the  arrival  time  of  the  So  mode  in  signal  So 
is  known,  the  non-zero  response  of  signal  So  outside  the  So  mode  is  mainly  attributed  to  Cso-  Therefore,  the  maximum 
amplitude  of  Cso  can  be  estimated  from  the  rest  of  signal  So  where  the  So  mode  does  not  exist.  The  maximum  absolute 
amplitude  of  cao  can  be  also  estimated  from  signal  Ao  in  a  similar  fashion.  Then,  the  threshold  value  for 
damage  classification  is  set  to  be  the  maximum  absolute  amplitude  of  cso  and  Cao-  Finally,  the  existence  of  damage  is 
determined  by  comparing  the  amplitudes  of  signals  MCl  and  MC2  with  respect  to  the  threshold  value.  In  this  way, 
a  crack  can  be  detected  if  the  crack  produces  signal  differences  between  signals  AB  and  CD  that  are  bigger 
than  the  initial  differences  due  to  PZT  imperfection  and  measurement  errors.  Note  again  that  no  predetermined 
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threshold  values  are  required  during  this  damage  classification  procedure  because  the  threshold  value  is 
obtained  from  instantaneously  measured  signals  AB,  AC,  BD,  and  CD.  The  applicability  of  the  proposed 
thresholding  technique  is  experimentally  investigated  in  the  following  subsection. 
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Figure  13.  Decomposition  of  signal  AB 
into  individual  Lamb  wave  mode  signals 
(When  a  notch  is  closer  to  PZT  A) 
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Figure  14.  Comparison  of  relative  phase  information 
among  signals  AB,  AC,  BD  and  CD  (When 
notch  is  closer  to  PZT  A) 

3.2.2.  Damage  diagnosis  using  experimentally  obtained  Lamb  wave  signals 

To  explain  the  proposed  damage  diagnosis  technique,  Lamb  wave  signals  measured  using  the  setup  in  3.1.2  were 
used.  First,  each  Lamb  mode  signal  is  estimated  from  signals  AB,  AC,  BD  and  CD  measured  from  an  unknown 
condition  of  the  specimen.  Figure  15  displays  individual  Lamb  mode  signals  decomposed  using  Equation  (2).  In 
Figure  15  (a),  the  estimated  signal  So  is  shown  as  a  superposition  of  the  exact  signal  So  and  the  error  signal,  Cso- 
Similarly,  signal  Ao  shown  in  Figure  15(b)  illustrates  that  the  error  signal  Cao  is  superimposed  to  the  exact  signal  Aq. 
Signals  MCl  and  MC2  are  also  drawn  in  Figure  15  (c)  and  Figure  15  (d).  Ideally,  signals  MCi  and  MC2  should  be 
zero  for  the  entire  length  of  the  signals  in  the  absence  of  damage.  However,  due  to  PZT  imperfections, 
additional  modes  appeared  in  signals  MCl  and  MC2. 

Once  the  individual  Lamb  mode  signals  are  decomposed,  the  threshold  for  damage  classification  is  set  to  be  the 
maximum  absolute  amplitude  of  error  components,  Cso  and  Cao.  in  signals  So  and  Aq.  The  estimated  signal  So  is 
composed  of  the  exact  signal  So  and  Csoas  shown  in  Figure  15.  Furthermore,  it  can  be  shown  that  the  So  mode  exists 
only  in  a  certain  portion  of  the  estimated  signal  So  while  Cso  is  present  for  the  entire  length  of  the  estimated  signal  So. 
Therefore,  the  maximum  absolute  amplitude  of  Cso  can  be  estimated  from  a  portion  of  the  estimated  signal  So  where 
the  So  mode  does  not  exist.  To  accomplish  this,  the  estimated  signal  So  is  divided  into  two  regions.  So  and  Ao  regions, 
as  shown  in  Figure  15.  Here,  the  boundary  between  two  regions  is  placed  exactly  in  the  middle  of  the  So  and  Ao 
modes.  Because  the  So  mode  always  exists  only  in  the  So  region,  the  response  in  the  Ao  region  is  attributed  only  to  Cso- 
Therefore,  the  maximum  absolute  amplitude  of  Cso  is  estimated  from  the  Ao  region  as  shown  in  Figure  15  (a).  In  a 
similar  manner,  the  maximum  absolute  amplitude  of  Cao  is  estimated  from  the  So  region  of  the  estimated  signals  Ao  as 
depicted  in  Figure  15  (b). 

From  the  estimated  maximum  absolute  amplitude  of  Cso  and  Cao,  the  upper  and  lower  limits  for  damage  diagnosis 
are  established  as  shown  in  Figure  16  (a).  It  is  assumed  that,  if  the  amplitudes  of  MCl  and  MC2  modes  exceed  one 
of  the  threshold  values,  the  mode  conversion  beyond  PZT  imperfections  or  measurement  noises  is  induced  by 
the  notch.  Note  that  these  threshold  values  for  damage  classification  are  strictly  estimated  from  signals  AB,  AC, 
BD  and  CD  instantaneously  measured  from  the  system’s  current  state.  That  is,  the  dependence  on  the  prior 
baseline  data  is  eliminated  during  the  damage  classification  process. 
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In  Figure  16  (a),  the  proposed  diagnosis  scheme  is  applied  to  the  signals  obtained  from  the  pristine 
condition  of  the  specimen.  Figure  16  (a)  illustrates  that  both  MCl  and  MC2  modes  in  the  So  and  Aq  regions  did  not 
exceed  the  threshold  values  determined  from  Cso  and  Cao-  Therefore,  it  can  be  concluded  that  these  signals  are 
obtained  from  the  undamaged  condition  of  the  test  article. 

Next,  the  proposed  scheme  was  applied  to  the  signals  measured  from  the  3  mm  notch  case.  Four 
decomposed  signals  using  Equation  (2)  are  drawn  in  Figure  15  along  with  the  signals  obtained  from  the  undamaged 
case  for  comparison.  Figure  15  (a)  and  Figure  15  (b)  showed  that  the  occurrence  of  damage  changed  signals 
So  and  Aq  very  little.  Furthermore,  it  is  revealed  that  the  amplitudes  and  shape  of  Cso  and  Cao  did  not  change 
much  because  they  were  not  related  to  mode  conversion.  On  the  other  hand,  the  appearance  of  mode  conversion  is 
clearly  observed  in  signals  MCl  and  MC2  [Figure  15  (c)  and  Figure  15  (d)].  After  the  threshold  values  are 
determined  from  the  maximum  absolute  amplitude  of  Cso  and  Cao,  it  is  tested  whether  the  MCl  rnode  in  the  Sq  region 
of  signal  MCl  and  the  MC2  mode  in  the  Aq  region  of  signal  MC2  exceed  the  thresholds.  Figure  16  (b)  shows  that 
both  MCl  and  MC2  modes  in  signals  MCl  and  MC2  exceeded  the  threshold  values  determined  from  the  damaged 
condition  of  the  structure.  Therefore,  these  signals  are  classified  to  be  damaged. 

In  this  section,  a  new  damage  classifier  that  does  not  rely  on  predefined  threshold  values  has  been  introduced. 
The  classifier  was  developed  based  on  the  decomposition  technique  introduced  in  3.2.1.  The  advantage  of  the 
decomposition  technique  is  that  only  converted  modes  can  be  extracted  from  measured  Lamb  wave  signals.  In  the 
following  parts  of  this  report,  this  decomposition  technique  plays  important  roles  for  damage  identification. 
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(c)  Decomposed  signal  MCl 
Figure  15.  Decomposition  of  individual  Lamb  mode  signals  Sq,  MCI,  MC2  and  Aq  from  measured  signals 
AB,  AC,  BD  and  CD  using  Equation  (2) 
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(a)  Comparison  of  the  amplitudes  of  estimated  MCi 
and  MC2  with  instantaneously  obtained  threshold 
values  (without  damage) 
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(b)  Comparison  of  the  amplitudes  of  estimated  MCl 
and  MC2  with  instantaneously  obtained  threshold 
values  (with  a  3  mm  notch) 
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Figure  16.  Autonomous  damage  identification  using  instantaneously  obtained  threshold  values  instead  of 
predetermined  decision  boundaries  (Max(abs(eso ,  Cao))  denotes  maximum  absolute  value  of  Cso  and  Cao- 


3.3.  Environmental  and  operational  variations 

In  this  project,  the  effect  of  environmental  and  operational  variations  on  the  proposed  technique  has  been  investigated. 
Especially,  experiments  were  conducted  in  the  controlled  laboratory  to  validate  the  robustness  of  the  proposed  NDT 
technique  to  varying  temperature  conditions. 

3.3.1.  Experimental  setup  for  varying  temperature  tests 

For  this  study,  an  aluminum  plate  of  455  mmx254  mmx3  mm  was  used  as  shown  in  Figure  17.  The  size 
of  the  plate  was  mainly  limited  by  the  available  space  of  the  temperature  chamber  used  later  for 
temperature  experiments.  The  Young’s  modulus  of  this  T6  aluminum  plate  was  310  MPa. 

Circular  PZTs  (6.35  mm  in  diameter  and  0.25  mm  in  thickness)  were  purchased  from  American  Piezo 
Ltd.  They  had  a  Curie  temperature  of  360^C  and  the  maximum  operational  temperature  of  180  ^C.  The  dss 
piezoelectric  charge  constant,  the  capacitance  value  and  the  Young’s  Modulus  were  4x  10'^^  m/V,  1.50  nF, 
6x  10^^  N/m^,  respectively.  PZTs  A  and  D  were  collocated  and  attached  on  the  other  side  of  the  plate,  and 
PZTs  B  and  C  were  mounted  in  a  similar  fashion.  The  PZTs  were  attached  so  that  their  poling  directions 
were  identical  to  the  configuration  shown  in  Figure  17.  (In  the  figure,  the  positive  poling  direction  of 
individual  PZT  is  shown  with  an  arrow.)  PZTs  A  and  B  (or  PZTs  C  and  D)  were  215  mm  apart  each  other. 

For  the  attachment  of  PZTs  to  the  specimen,  M-bond  200cyanoacrylate  adhesive  from  Vishay 
measurements  group  was  used.  Based  on  the  manufacturer’s  specification,  this  M-bond  adhesives  could  be 
used  for  a  one-cycle  proof  test  over  90  or  -185  ^C.  But,  the  normal  operating  temperature  range  was  -30 
to  65  ^C.  To  secure  the  wires  to  the  specimen,  high  temperature  Teflon  tapes  from  3M  Corporation  (3M 
Scotch  Brand  5412)  were  used.  In  addition,  high  temperature  wires  were  used. 


Figure  17.  Dimension  and  configuration  of  the  aluminum  test  specimen  with  a  uniform  thickness: 
The  arrows  denote  the  positive  polarization  directions  of  each  PZTs. 


AOARD07  FINAL  REPORT 


The  data  acquisition  system  used  in  this  study  is  shown  in  Figure  18.  The  data  acquisition  system  was 
composed  of  a  laptop  computer,  an  arbitrary  waveform  generator  (AWG),  a  high-speed  signal  digitizer  (DIG),  a  low 
noise  preamplifier  (LNP)  and  a  multiplexer.  Using  the  14-bit  AWG,  a  toneburst  signal  with  a  10  peak-to-peak 
voltage  and  a  driving  frequency  of  280kHz  was  generated  and  applied.  First,  PZT  A  shown  in  Figure  17  was 
excited  by  this  input  waveform.  Then,  PZT  A  generated  elastic  waves  and  the  response  was  measured  at  PZT  B. 
When  the  waves  arrived  at  PZT  B,  the  voltage  output  from  PZT  B  was  filtered  and  amplified  by  the  LNP  with  a  gain 
of  ten  and  measured  by  the  DIG.  The  sampling  rate  and  resolution  of  the  DIG  were  20  MS/sec  and  16  bits, 
respectively.  In  order  to  improve  the  signal-to-noise  ratio,  the  forwarding  signals  were  measured  ten  times  and 
averaged  in  the  time  domain.  After  the  forwarding  signals  from  PZT  A  to  PZT  B  (signal  AB)  were  measured,  the 
same  process  was  repeated  by  measuring  signals  AC,  DB  and  DC.  The  entire  experimental  process  without  averaging 
took  less  than  1  minute.  By  increasing  the  number  of  channels,  it  is  expected  that  the  data  collection  time  can  be 
further  shortened.  The  switching  of  the  exciting  and  sensing  PZTs  was  accomplished  using  the  multiplexer  and 
controlled  by  the  Lab  VIEW  software  installed  on  the  laptop  computer.  Initial  data  was  collected  without  any  crack  in 
the  specimen,  and  then  a  crack  was  introduced  between  PZTs  A  and  B,  54  mm  from  PZT  A  as  shown  in  Figure  19. 
The  length  and  width  of  the  crack  was  3 1  mm  and  1  mm,  respectively,  and  the  crack  depth  was  increased  from  0  mm 
to  0.5  mm,  1.0  mm  and  2.0  mm. 

As  shown  in  Figure  17,  two  thermocouples  T1  and  T2  are  instrumented  at  the  middle  and  quadrant 
points  of  the  specimen.  The  thermocouples  were  connected  to  BK  Precision  725  Humidity  and  Temperature  Meter. 
As  part  of  temperature  experiments  shown  in  Figure  20(a),  an  infrared  heater  from  Protherm  Heater  (model  #12014) 
was  used  to  simulate  partial  head-up  of  the  specimen.  For  the  experimental  results  provided  in  the  following 
subsection,  the  reading  from  thermocouple  T2  was  reported. 

Another  temperature  experiment  was  conducted  using  a  MicroClimate  temperature/humidity  chamber.  The 
chamber  is  designed  to  operate  in  a  temperature  range  of  -73''C  to  lOO^'C  and  a  maximum  relative  humidity 
of  95%.  In  this  study,  only  the  temperature  value  was  controlled. 
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(a)  A  inferred  heater  (b)  A  temperature  chamber 

Figure  20.  Two  experimental  setups  for  temperature  variations 
3.3.2.  Crack  detection  under  varying  temperatures 

During  the  experiments,  multiple  Lamb  wave  signals  are  measured  from  the  aluminum  plate  in  both  intact 
and  damaged  conditions  under  varying  temperatures.  Since  the  size  of  the  plate  was  not  large  enough  to  avoid 
reflective  Lamb  wave  modes  from  the  boundary,  the  damage  classifier  proposed  in  section  3.2  could  not  be  used 
in  this  specific  case.  Note  that  the  error  signals,  which  the  former  classifier  utilized,  were  superimposed  by  the 
reflections  in  this  case.  Therefore,  a  new  thresholding  technique  (Technique  2)  was  proposed  in  this  subsection. 

Technique  2  starts  by  decomposing  measured  signals  AB,  AC,  DB  and  DC  into  signals  Sq,  MCi,  MC2 
and  Aq  again.  Then,  the  standard  deviation  of  signal  MCl,  Omci  is  computed.  If  there  is  no  defect,  the  variance  of 
signal  MCl  is  mainly  attributed  to  the  initial  errors.  Assuming  that  the  initial  errors  have  a  normal  distribution  of  zero 
mean  and  Omci  standard  deviation,  99.7%  of  the  data  points  in  signal  MCl  should  be  within  the  range  of  -3aMci 
to  3aMci-  Once  actual  mode  conversion  appears  in  signal  MCl,  the  maximum  absolute  amplitude  of  the 
actual  MCl  mode  should  go  outside  this  range.  Therefore,  the  existence  of  a  crack  can  be  identified  when  max(MCl) 
becomes  larger  than  Omci-  In  a  similar  fashion,  the  threshold  for  signal  MC2  can  be  set  to  be  aMC2- 


Ifmax([MC^  |)  mid  max(|IV[C2 1)  >  3c7y^^2^  a  crack  exists 


(3) 


The  next  step  was  to  examine  if  the  proposed  NDT  technique  and  damage  classifier  are  robust  even  under  changing 
temperature  conditions.  The  variations  of  measured  time  signals  and  decoupled  Lamb  modes  with  respect  to 
temperature  were  illustrated  in  Figure  21  using  representative  signals.  These  signals  were  obtained  when  the 
specimen  with  2mm  deep  crack  was  placed  inside  the  temperature  chamber.  The  amplitude  increase  and  phase  delay 
of  the  signals  were  observed  as  the  temperature  increased  from  -30  ‘'C  to  70  ""C.  It  should  be  noted  that  because 
the  temperature  variation  itself  could  cause  significant  changes  in  the  measured  signal’s  amplitude  and  phase, 
it  would  be  challenging  to  perform  pattern  comparison  with  the  baseline  signals  for  the  purpose  of  damage  diagnosis. 

All  damage  diagnosis  results  from  varying  temperature  conditions  are  summarized  in  Table  2.  In  the  first  4  rows 
(cases  1-4),  false-positive  studies  were  conducted  using  the  intact  specimen.  For  all  cases  examined,  there  were  no 
indications  of  false  alarms.  For  cases  1  to  4,  the  temperature  was  controlled  using  the  infrared  heater  from  22  ""C  to  55 
‘'C.  Note  that  the  temperature  reported  here  was  obtained  from  thermocouple  T2  shown  in  Figure  17,  and 
thermocouple  T1  was  used  as  the  set  point  for  the  temperature  controller.  Because  the  infrared  device  locally 
heated  up  the  middle  portion  of  the  specimen,  a  large  temperature  gradient  between  the  center  and  the  edge 
of  the  specimen  was  observed.  For  instance,  the  reading  form  T2  was  close  to  100  '’C  for  case  4  in  Table  2,  while  T1 
was  about  50  ""C.  Note  that  a  broader  spectrum  of  temperature  variations  was  not  examined  using  the  intact 
specimen  to  avoid  any  potential  damage  to  the  PZT  transducers  before  performing  subsequent  damage  cases.  For 
cases  5  to  16  in  Table  2,  the  specimen  with  1  mm  crack  was  subjected  to  temperature  increase  from  5''C  to  bO^'C  using 
either  the  infrared  heater  or  the  temperature  chamber.  For  all  cases  investigated,  the  maximum  absolute  amplitudes  of 
the  MCl  and  MC2  modes  exceeded  the  threshold  values  obtained  from  Technique  2.  Furthermore,  the  numbers  of 
outliers  for  each  threshold  values  remained  reasonably  consistent.  For  the  remaining  cases,  (cases  17  -  27),  the 
specimen  with  2  mm-deep  crack  was  subjected  a  broad  temperature  variation  of -30®C  to  70  using  the  temperature 
chamber.  Compared  to  cases  5  to  17,  larger  numbers  of  outliers  were  observed  for  cases  17  to  27. 

It  was  concluded  that  a  crack  deeper  than  1  mm  was  detectable  when  the  3  mm  specimen  was  investigated  in  the 
temperature  range  of  -30  '’C  to  70  ""C  as  long  as  there  was  no  degradation  of  the  PZTs’  piezoelectricity  and  bonding 
conditions.  The  applicability  of  the  proposed  technique  was  demonstrated  under  the  limited  conditions:  The 
consistent  bonding  of  the  PZTs  and  the  precise  alignment  of  the  collocated  PZTs  were  critical  to  the  success  of 
the  proposed  technique.  In  addition,  the  current  technique  is  only  applicable  to  structures  with  uniform 
thickness.  As  long  as  mode  conversions  reflected  or  refracted  from  the  crack  are  measurable  by  one  of  the 
PZTs,  the  mode  conversion  due  to  the  crack  formation  should  be  detected  in  theory.  However,  the  effects  of 
the  crack  location  and  orientation  were  not  investigated  in  this  study.  Ongoing  research  is  underway  to  address  these 
issues. 
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(a)  Measured  signal  AB 


0.02  0.04 

Time  (ms) 

(b)  Measured  signal  BC 


0.06 


(c)  Decoupled  signal  So 


(d)  Decoupled  signal  Aq 
■-30°C  0°C  — ■ 


70°C 


(e)  Decoupled  signal  MCi  (Sq/Ao) 


0.02  0.04 

Time  (ms) 

(f)  Decoupled  signal  MC2  (Aq/So) 


0.06 


Figure21.  Variation  of  measured  time  signals  and  decomposed  modes  with  repect  to  temperature  ( with  2 
mm  deep  crack :  -30°C  0°C  22°C  70°C ) 


Table  2:  Number  of  outliers  outside  the  threshold  values  (under  changing  temperatures) 


Crack 

depth 

# 

T2  (°C)* 

Technique  2 

Temp. 

Control 

Date 

(mm/dd/yy) 

Sq/Aq 

Aq/Sq 

0  mm 

1 

22 

0 

0 

Infrared 

heater 

08/01/07 

2 

39 

0 

0 

08/01/07 

3 

40 

0 

0 

08/01/07 

4 

55 

0 

0 

08/01/07 

1  mm 

5 

23 

47 

31 

Infrared 

heater 

08/07/07 

6 

30 

44 

39 

08/07/07 

7 

40 

44 

39 

08/07/07 

8 

51 

44 

39 

08/07/07 

9 

05 

43 

32 

Temp. 

chamber 

08/07/07 

10 

10 

45 

28 

08/08/07 

11 

15 

44 

31 

08/08/07 

12 

21 

45 

31 

08/09/07 
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13 

23 

44 

31 

08/08/07 

14 

41 

45 

37 

08/09/07 

15 

50 

42 

36 

08/09/07 

16 

60 

42 

40 

08/09/07 

17 

-30 

43 

41 

08/13/07 

18 

-20 

46 

43 

08/13/07 

19 

-10 

45 

41 

08/13/07 

20 

-5 

49 

49 

08/13/07 

21 

0 

47 

39 

Temp. 

chamber 

08/13/07 

2  mm 

22 

11 

48 

43 

08/13/07 

23 

17 

47 

41 

08/13/07 

24 

22 

46 

46 

08/13/07 

25 

50 

47 

47 

08/13/07 

26 

60 

51 

55 

08/13/07 

27 

70 

50 

51 

08/13/07 

*The  temperature  reported  here  was  obtained  from  thermocouple  T2  shown  in  Figure  17. 


3.4.  Extension  to  more  complex  structures 


So  far,  the  proposed  technique  has  been  tested  for  simple  plates  although  many  civil  infrastmctures  have  more 
complex  stmctural  features  such  as  holes  and  stiffeners.  In  this  project,  not  only  the  effect  of  through-the- 
thickness  holes  but  also  stiffeners  on  the  proposed  technique  is  investigated. 


3.4.1.  Understanding  the  effect  of  through-the-thickness  holes 

An  array  of  holes  creates  multiple  reflections  and  refractions  of  Lamb  waves,  and  these  multiple  reflections  and 
refractions  may  cause  difficulties  in  analyzing  responses  of  Lamb  waves  [5].  Because  a  symmetric 
configuration  of  waveguide  thickness  variation  is  preserved  even  in  the  existence  of  through-the-thickness  holes, 
the  holes  do  not  produce  mode  conversion  [2].  Therefore,  it  is  expected  that  those  holes  do  not  affect  the 
performance  of  the  proposed  technique.  This  theoretical  expectation  was  validated  from  numerical  simulations  and 
experiments. 

Since  the  through-the-thickness  hole  could  not  be  modeled  under  plain  strain  assumption.  Lamb  wave 
propagation  in  a  three-dimensional  aluminum  plate  was  simulated  using  PZFlex  1.4.  The  configuration  of  specimen 
including  a  hole,  a  notch  and  PZT  wafer  transducers  are  shown  in  Figure  22. 


Figure  22.  Configuration  of  the  specimen  with  collocated  PZTs,  a  notch  and  a  through-the-thickness  hole  for 
numerical  simulation  (all  dimension  in  mm).  The  arrows  indicate  the  positive  poling  direction  of  PZTs 

In  Figure  23,  simulation  results  are  illustrated.  As  expected,  even  in  the  presence  of  the  through-the-thickness  hole, 
simulated  signals  AB  and  DC  (or  signals  AC  and  DB)  were  identical  because  the  hole  was  not  the  source  of  mode 
conversion  [Figure  23  (a)  and  Figure  23  (b)].  Signals  AB  and  DC  (or  signals  AC  and  DB)  became  different  after  the 
notch  was  formed  [Figure  23  (c)  and  Figure  23  (d)].  These  results  are  similar  to  the  simple  plate  case. 
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X  10“^  X  ^0‘^ 


(a)  Signals  AB  and  DC  with  only  a  hole 


(b)  Signals  AB  and  DC  with  a  hole  and  a  notch 


X  10“^  X  10"^ 


(c)  Signals  AC  and  DB  with  only  a  hole 


(d)  Signals  AC  and  DB  with  a  hole  and  a  notch 


Figure  23.  Comparison  of  simulated  signals  AB  and  DC  (or  signals  AC  and  DB) 

After  So,  Aq,  MCl  and  MC2  were  decoupled  using  Equation  (2),  it  was  shown  clearly  that  the  through-the- 
thickness  hole  did  not  produce  converted  modes  [Figure  24  (c)  and  Figure  24  (d)].  Also,  the  patterns  of  the  So, 
Aq,  MCl,  and  MC2  modes  did  not  change  much  even  in  the  presence  of  the  hole. 


(a)  Decomposed  So  modes 


(b)  Decomposed  Ao  modes 


(c)  Decomposed  MCi  modes 


(d)  Decomposed  MC2  modes 
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- With  only  a  hole  - With  a  hole  and  a  notch - With  only  a  notch 

Figure  24.  Decomposed  signals  So,  MCi,  MC2  and  Aq  from  simulated  Lamb  wave  signals  using  Equation  (2), 
where  signal  MCi  and  MC2  represent  the  first  and  the  second  mode  conversion,  respectively. 

After  the  numerical  simulation,  the  effect  of  the  through-the-thickness  holes  has  been  investigated 
experimentally.  The  aluminum  plate  used  in  this  study  has  identical  dimensions  with  the  one  in  Figure  22.  Except  for 
the  input  frequency  (180kHz)  and  the  gain  of  LNP  (50),  other  parameters  were  selected  to  be  identical  to  the  values  in 
section  3.1. 

First,  one  hole  and  the  notch  were  introduced  in  series  and  the  responses  from  that  specimen  were  measured 
[Figure  25].  Next,  two  additional  holes  were  formed  and  similar  procedure  was  repeated.  The  three  holes  were  50 
mm  apart  to  each  other.  The  diameter  of  the  holes  and  the  dimensions  of  the  notch  are  shown  in  Figure  25  (b). 


(a)  The  specimen  with  collocated  PZTs,  a  notch  and  (b)  A  notch  was  produced  first  and  then  a  through-the- 
through-the-thickness  holes  for  experiments.  thickness  hole  was  introduced 

Figure  25.  The  specimen  for  experiments  to  consider  the  effect  of  the  through-the-thickness  holes 

Figure  26(a)  and  Figure  26(b)  clearly  showed  that  the  shapes  and  amplitudes  of  the  unconverted  modes  (Sq/  Sq 
and  Aq/Aq  modes)  were  not  affected  much  by  either  the  notch  or  the  holes.  Additionally,  converted  Lamb 
wave  modes  due  to  the  notch  appeared  in  signals  MCi  and  MC2  regardless  to  the  number  of  holes.  Also  note  that,  the 
damage  classifier  proposed  in  section  3.2  can  be  applied  to  these  decomposed  signals  because  the  effect  of  the 
through-the-thickness  holes  on  the  proposed  technique  is  minimal. 


(a)  Comparison  of  signals  Sq 


(b)  Comparison  of  signals  Aq 
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(c)  Comparison  of  signals  MCi  (d)  Comparison  of  signals  MC2 

- With  only  a  hole  - With  a  hole  and  a  notch - With  3  holes  and  a  notch 

Figure  26.  Comparison  of  decomposed  signals  among  the  condition  of  the  plate  with  only  a  hole,  with  a  hole 
and  a  notch  and  with  3  holes  and  a  notch 

3.4.2.  Understanding  the  effect  of  stiffeners 

In  this  subsection,  the  proposed  damage  detection  technique  is  extended  to  the  plate  including  a  stiffener. 
Since  the  configuration  of  waveguide  become  no  longer  symmetrical  when  the  stiffener  is  introduced  to  the  wave 
path,  mode  conversion  can  occur  even  though  damage  is  not  formed.  Mode  conversion  due  to  the  stiffener  may 
cause  difficulties  in  analyzing  the  Lamb  wave  signals  measured  from  the  structures. 

In  Figure  27,  some  of  the  converted  Lamb  wave  modes  due  to  the  stiffener  and  notch  are  illustrated. 
In  the  plate  with  a  single  stiffener,  there  exist  two  possible  wave  paths:  (1)  the  wave  propagation  only  inside  the 
plate  [Figure  27  (a)],  and  (2)  the  wave  path  through  the  stiffener  [Figure  27  (b)].  In  the  first  path,  a 
single  mode  is  separated  into  four  individual  modes  because  there  are  two  sources  of  mode  conversion, 
the  stiffer  and  the  notch.  (Underline  denotes  that  the  mode  is  produced  by  the  notch.)  Therefore,  the  So  and  Aq 
modes  produce  total  8  modes  in  the  first  path.  In  the  second  path,  three  sources  of  mode  conversion  exist.  The  first 
and  the  second  mode  conversion  are  occurring  when  the  single  mode  is  entering  and  propagating  out  of  the  stiffer. 
The  third  source  of  mode  conversion  is  the  notch.  At  the  end,  a  single  Lamb  wave  mode  is  separated  into  8 
individual  modes  in  the  second  path.  Finally,  total  24  Lamb  wave  modes  are  measured  at  the  receiving  transducer 
in  the  presence  of  the  stiffener  and  the  notch. 
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(a)  Converted  and  transmitted  Lamb  modes  due  to  (b)  Converted  and  reflected  Lamb  modes  due  to  a 

a  stiffener  and  a  notch  on  the  direct  path  stiffener  and  a  notch  on  the  path  including  a 

stiffener 

Figure  27.  A  scheme  of  Lamb  wave  propagation  on  the  plate  including  a  stiffener 

The  PI  is  now  preparing  for  a  journal  paper  that  includes  detailed  analysis  of  each  Lamb  wave  modes  in  the  plate 
with  one  stiffener  case.  The  main  idea  of  the  journal  paper  is  introduced  briefly  in  this  report.  From  measured  Lamb 
wave  signals,  signals  MCi  and  MC2  can  be  decomposed.  In  the  absence  of  damage,  these  two  signals  are  identical 
except  that  signal  MCl  is  shifted  compared  to  signal  MC2.  This  shifting  in  two  signals  is  caused  when  the  stiffener  is 
not  exactly  in  the  middle  between  PZTs  A  and  B  (or  PZTs  D  and  C).  Also,  using  information  about  group  velocities 
of  the  So  and  Aq  modes  and  the  geometry  of  the  structure,  the  amount  of  shifting  can  be  estimated.  Therefore,  the 
mode  conversion  produced  by  only  stiffener  can  be  compensated  by  aligning  two  signals. 
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On  the  other  hand,  signals  MCi  and  MC2  are  different  even  after  the  alignment  in  the  presence  of  the  notch.  Figure 
28  illustrates  each  mode  in  two  signals.  For  example,  even  thought  the  SSA  mode  in  shifted  signal  MCl  and 
the  ASS  mode  in  signal  MC2  have  same  arrival  times,  they  have  different  amplitudes  and  shapes.  When  signal  MC2 
is  subtracted  from  the  shifted  MCl,  multiple  Lamb  wave  modes  will  remain.  This  idea  was  fully  validated  by 
numerical  simulations  as  well  as  experiments. 
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Figure  28.  Extracted  MCi  and  MC2  modes  of  measured  signals  from  the  damaged  specimen 

In  Figure  29,  the  dimensions  of  the  aluminum  plate  with  one  stiffener  used  in  this  experiment  are  shown.  In  the 
simulation,  the  plate  Figure  29  was  modeled  as  a  2-D  stmcture  under  plain  strain  assumption  for  brevity.  COMSOL 
3.4  Multiphysics  software  was  used  to  simulate  Lamb  wave  propagation  in  a  two  dimensional  aluminum  plate 
with  one  stiffener.  The  length  of  the  plate  was  70  cm,  and  its  thickness  was  6  mm.  The  depth  and  the  thickness  of  the 
stiffener  were  90  mm,  and  6  mm,  respectively.  Four  identical  PZTs  with  10  mm  length  and  0.5  mm  thickness  were 
attached  to  the  plate  model  as  shown  in  Figure  29.  The  parameters  in  Table  1  were  used  in  this  simulation. 

Signals  MCl  and  MC2  decomposed  from  the  Lamb  wave  signals  are  shown  in  Figure  30.  Figure  30  (a) 
and  Figure  30  (c)  illustrate  that  the  mode  conversion  due  to  only  the  stiffener  can  be  compensated  by  shifting  signal 
MCl  in  the  absence  of  damage  [Figure  30  (e)].  Once  a  notch  of  2  mm  depth  and  1mm  width  was  introduced  right 
next  to  the  stiffener  as  is  shown  in  Figure  29,  shifted  signal  MCl  became  different  from  signal  MC2  as  a  result 
of  the  mode  conversion  induced  by  the  notch  [Figure  30  (d)  and  Figure  30  (f)]. 


700 


Figure  29.  Configuration  of  the  specimen  with  collocated  PZTs,  a  notch  and  a  stiffener  for  numerical 
simulation  (all  dimension  in  mm).  The  arrows  indicate  the  positive  poling  direction  of  PZTs 


AOARD07  FINAL  REPORT 


(a)  Comparison  of  signals  MCi  and  MC2  (with  only  a 
stiffener) 


(c)  Comparison  of  signals  MCi  and  MC2  (with  a 
stiffener  and  a  notch) 
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(e)  Difference  of  MC2  from  shifted  MCi  (with  only  a 
stiffener) 


(b)  Consideration  of  the  difference  of  arrival  time 
between  signals  MCi  and  MC2  (with  only  a 
stiffener) 


(d)  Consideration  of  the  difference  of  arrival  time 
between  signals  MCi  and  MC2  (with  a  stiffener  and 
a  notch) 


(f)  Difference  of  MC2  from  shifted  MCi  (with  a  stiffener 
and  a  notch) 


—  Shifted  MC1  -  MC2 


Figure  30.  Comparison  of  individual  Lamb  modes  and  compensation  of  the  difference  in  the  arrival  times 
between  MCi  and  MC2  using  simulated  Lamb  wave  signals 


After  the  numerical  simulation,  the  effect  of  the  stiffener  has  been  investigated  experimentally.  The  aluminum  plate 
used  in  this  study  has  identical  dimensions  with  the  one  in  Figure  29.  Except  for  the  gain  of  LNP  (50),  other  parameters 
were  selected  to  be  identical  to  the  values  in  section  3.1. 
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(a)  Configuration  of  a  plate  including  a  stiffener  and  PZTs 


(b)  The  dimension  of  a  stiffener 
Figure  31.  Configuration  of  the  specimen 
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(c)  The  location  and  dimension  of  a  notch 


Signals  MCl  and  MC2  decomposed  from  the  measured  Lamb  wave  signals  are  shown  in  Figure  32.  In  the 
simulation  results,  signals  MCl  and  MC2  were  identical  after  alignment  when  the  notch  was  absent.  On  the  other 
hand,  Figure  32(c)  showed  that  the  two  signals  had  different  amplitudes.  It  was  suspected  that  the  difference  in 
amplitudes  was  caused  by  the  imperfection  in  PZT  bonding  conditions  and  different  welding  conditions  on  both  sides 
of  the  stiffener.  However,  signals  MCl  and  MC2  were  almost  identical  in  shape.  Figure  32(d)  showed  that  signals 
MCl  and  MC2  were  different  in  both  amplitude  and  shape  due  to  the  presence  of  the  notch.  As  a  part  of  the  project, 
the  damage  classifier  for  this  plate  with  one  stiffener  case  is  now  being  developed  and  tested. 


(a)  Comparison  of  signals  MCi  and  MC2  (with  only  a  (b)  Comparison  of  signals  MCi  and  MC2  (with  a 
stiffener)  stiffener  and  a  notch) 
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(c)  Consideration  of  the  difference  of  arrival  time 
between  signals  MCi  and  MC2  (with  only  a 
stiffener) 


(d)  Consideration  of  the  difference  of  arrival  time 
between  signals  MCi  and  MC2  (with  a  stiffener  and 
a  notch) 


Figure  32.  Comparison  of  individual  Lamb  modes  and  compensation  of  the  difference  in  the  arrival  times 
between  MCi  and  MC2  using  measured  Lamb  wave  signals 


4.  Summary  and  Conclusion 


In  this  project,  the  applicability  of  the  proposed  NDT  technique  to  crack  detection  has  been  fully  validated  throughout 
numerical  as  well  as  experimental  studies.  The  PI  listed  the  specific  objectives  of  the  project  below  and  elaborated 
on  the  findings  in  each  objective. 

Extraction  of  Reference-Free  Features:  In  the  project,  the  idea  of  using  PZT  polarization  characteristics  for 
extracting  damage  features  has  been  tested  numerically  and  experimentally.  The  phase  information  of  the 
measured  Lamb  wave  signals  in  a  time  domain  depends  on  whether  the  sensing  PZT  is  attached  to  the  top  or  the 
bottom  surface  of  the  plate.  Furthermore,  the  phases  of  the  generated  Lamb  waves  change  in  accordance  with  the 
poling  direction  of  the  exciting  PZT  transducers.  Therefore,  by  controlling  the  poling  directions  of  the  exciting 
and  sensing  PZT  transducers,  it  is  possible  to  collect  different  Lamb  wave  signals  in  a  wave  path.  In  this  project,  the 
use  of  four  collocated  PZT  transducers  was  proposed  to  measure  multiple  Lamb  wave  signals.  When  propagating 
Lamb  waves  experienced  mode  conversion  in  the  presence  of  damage,  only  the  converted  modes  due  to  the 
damage  could  be  successfully  extracted  from  the  Lamb  wave  signals  measured  using  the  PZT  transducers. 
Instantaneous  Damage  Diagnosis:  Once  damage-sensitive  features  were  identified,  a  theoretical  framework  for 
statistical  damage  diagnosis  was  developed  to  fully  automate  the  decision-making  procedure.  The  proposed  damage 
diagnosis  aimed  to  differentiate  signals  due  to  damage  from  error  signals  produced  by  variations  in  PZT  size, 
alignment,  and  bonding  condition.  In  this  project,  it  has  been  shown  that  the  measured  Lamb  wave  signals  could  be 
decomposed  so  that  converted  modes  due  to  damage  could  be  extracted.  The  decision  boundary  was  chosen  from  the 
error  signals  extracted  only  from  the  measured  Lamb  wave  signals.  By  comparing  the  amplitudes  of  the  converted 
modes  with  those  of  error  signals,  automated  decision  was  made.  Note  that  only  instantaneously  measured  Lamb 
wave  signals  were  used  to  identify  the  current  condition  of  the  structure. 

Environmental  and  Operational  Variations:  Varying  temperature  tests  have  been  conducted  in  the  controlled 
laboratory.  Four  PZT  transducer  wafers  were  attached  to  an  aluminum  plate,  and  converted  modes  due  to  damage 
were  extracted  from  measured  Lamb  wave  signals.  As  long  as  the  PZTs’  piezoelectricity  and  bonding  conditions 
were  preserved,  the  crack  damage  was  detectable  in  the  temperature  range  of  -30  ‘'C  to  70  ‘'C. 

Extension  to  More  Complex  Structures:  Once  the  basic  concept  of  the  reference-free  NDT  technique  was 
established,  it  has  been  extended  for  more  complex  structures  such  as  a  plate  with  holes  and  a  plate  with  one  stiffener. 
Since  the  through-the-thickness  holes  were  not  the  sources  of  mode  conversion,  the  proposed  technique  could  be 
successfully  applied  to  the  plate  with  multiple  holes.  As  for  the  plate  with  one  stiffener,  it  has  been  shown  that  mode 
conversion  regardless  to  the  damage  could  be  compensated  so  that  only  the  converted  modes  due  to  the  damage  could 
be  identified  using  the  proposed  technique. 
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